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Back-to-back hadron pair yields in d+Au and p+p collisions at ^/snn = 200 GeV were measured 
with the PHENIX detector at the Relativistic Heavy Ion Collider. Rapidity separated hadron pairs 
were detected with the trigger hadron at pseudorapidity \n\ < 0.35 and the associated hadron at 
forward rapidity (deuteron direction, 3.0 < n < 3.8). Pairs were also detected with both hadrons 
measured at forward rapidity; in this case the yield of back-to-back hadron pairs in d+Au collisions 
with small impact parameters is observed to be suppressed by a factor of 10 relative to p+p collisions. 
The kinematics of these pairs is expected to probe partons in the Au nucleus with a low fraction 
x of the nucleon momenta, where the gluon densities rise sharply. The observed suppression as a 
function of nuclear thickness, pr, and n points to cold nuclear matter effects arising at high parton 
densities. 

PACS numbers: 25.75.Dw 



Nuclear effects on quark and gluon densities of nu- 
cleons bound in nuclei can be studied in collisions of 
deuteron and gold nuclei at the Relativistic Heavy Ion 
Collider. At ^smn = 200 GeV, the hadron yields in 
the forward-rapidity (deuteron-going) direction were ob- 
served to be suppressed for d+Au collisions relative to 
p+p collisions [lH3j. However, the mechanism for the 
suppression was not firmly established and may indicate 
novel QCD effects in nuclei. Competing theoretical ap- 
proaches include initial state energy loss 0, [![, parton 
recombination 6], shadowing effects [7, 8], and gluon sat- 
uration Q. 

Back-to-back dijet yields were proposed as an addi- 
tional observable to distinguish better between compet- 
ing mechanisms. The color glass condensate (CGC) 
framework predicts that quarks and gluons scatter- 
ing at forward angles (large rapidity) will interact coher- 
ently off gluons at low x in the gold nucleus. As a result, 
the rate of observed recoiling jets is expected to be sup- 
pressed in d+Au collisions compared to p+p, and angular 
broadenin g of the back-to-back correlation of jets is pre- 
dicted Jn|jl3. Dihadro n correlation measurements were 
used 13 , 14| successfully in d+Au collisions to select di- 
jet production based on the back-to-back peak at A4>—tt 
between trigger hadrons and their associated hadrons. 
Dihadron correlation measurements with varying kine- 
matic constraints (transverse momentum px and rapid- 
ity) probe different x ranges in the nucleus. In particular, 
measurements at forward rapidity are thought to probe 
small x values in the Au nucleus. 

In this Letter, we report results on the suppression 
in d+Au relative to p+p collisions of inclusive 7r°'s and 
back-to-back cluster-7r° pairs in the forward-rapidity re- 
gion, and for back-to-back ir°-ir° or hadron-7r° pairs 
separated in rapidity. The data were obtained from 
p+p and d+Au runs in 2008 with the PHENIX detec- 
tor and include a new electromagnetic calorimeter, the 
muon piston calorimeter (MPC), with an acceptance of 



3.0 < T] < 3.8 in pseudorapidity and < < 2n. The 
clusters are reconstructed from the energy deposit of pho- 
tons in individual MPC towers. The MPC comprises 220 
PbWC>4 towers of 20.2A+, depth, with lateral dimensions 
of 2.2x2.2 cm 2 , and is located 220 cm along the beam 
axis from the nominal interaction point. 

The d+Au sample is separated into four centrality 
classes - 0-20% (most central), 20-40%, 40-60%, and 
60-88% (most peripheral) - based on charge deposited 
in the backward (gold direction) beam-beam counter 
(3.0 < —r\ < 3.9). We determine the average number 
of binary collisions (iV co ii) from a Glauber model Q and 
a simulation of the beam-beam counter; (N co u) values 
are 15.1 ± 1.0, 10.2 ± 0.70, 6.6 ± 0.44, and 3.2 ± 0.19, 
respectively. 

The charged hadron (/i*) and ir° analysis in the midra- 
pidity region \r]\ < 0.35 is identical to that for previous 
measurements by PHENIX [TBI , [lij . For the analysis in 
the forward-rapidity region, a fiducial cut is applied en- 
suring that clusters are fully reconstructed within the 
MPC acceptance. Photon candidates are identified in 
the MPC by comparing cluster candidates to the ex- 
pected shower profile for photons. The shower profile for 



photons is determined from the PHENIX GEANT3 [17 1 
based detector simulation, PISA, which was tuned to re- 
produce MPC test beam data. For the pair analyses re- 
ported here, the associated particles in the MPC are 7r°'s, 
which are identified by reconstructing the mass of their 
decay photon pairs. The 7r° yield is obtained after sub- 
traction of the combinatoric background of uncorrelated 
photon pairs. The shape of the background was deter- 
mined from p+p pythia 6.4 18] and d+Au HIJING [ljj[ 
events that are subsequently processed through PISA. 
The pT-dependent systematic uncertainty on the asso- 
ciated 7T° yield extraction procedure is estimated to be 
2-5% for p+p and d+Au. 

The closeness of the MPC to the collision vertex and 
the high energy of particles emitted in the forward direc- 
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FIG. 1: (color online). Pedestal-subtracted 7r°-7r° per-trigger 
correlation functions for, as indicated, p+p, d+Au periph- 
eral (60-88% centrality) and d+Au central (0-20% central- 
ity) collisions at ^/snn = 200 GeV; the associated 7r°'s 
of pr = 0.5-0.75 GeV/c are measured at forward rapidity 
(3.0 < r) < 3.8) and the triggered 7r°'s are measured at 
midrapidity (\rj\ < 0.35) for the indicated pr ranges. The 
subtracted pedestal values, bo, are also indicated. 
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FIG. 2: (color online). Pedestal-subtracted cluster- tv° per- 
trigger correlation functions measured at forward rapidity 
(3.0 < r\ < 3.8) for, as indicated, p+p, d+Au peripheral 
(60-88% centrality) and d+Au central (0-20% centrality) col- 
lisions at ^snn = 200 GeV; the correlation functions are for 
associated 7r°'s of pr ~ 0.5-0.75 GeV/c and trigger clusters 
over the indicated pt ranges. Systematic uncertainties of up 
to 30% on the near side (|A0| < 0.5) are not shown. The 
subtracted pedestal values, bo, are also indicated. 



tion make it difficult to reconstruct photon pairs from 7r° 
decays at high pt- For example at Pt = 1 GeV/c, ap- 
proximately 30% of the photon cluster pairs are merged 
and cannot be reconstructed separately in the MPC. To 
extend the pt range and the pair yield, single electromag- 
netic clusters are used as trigger particles to construct 
cluster- 7T° dihadron pairs in the MPC. These trigger clus- 
ters are treated assuming that they are all 7r°'s. How- 
ever, pythia studies indicate that > 80% of these trig- 
ger clusters are from 7r° 's with the rest being dominantly 
single photons from asymmetric decays of r\ mesons or 
direct photons; thus, according to these studies a rela- 
tively small contamination remains. The cluster energy 
was corrected to the true 7r° energy to account for the 
merging effects of the two photons from 7r° decay. These 
corrections were determined by embedding Monte Carlo 
generated 7r°'s into real data, as well as from pythia 



tuned to match the data. 

Figure [T] shows the azimuthal angle correlations be- 
tween midrapidity and forward-rapidity tt° pairs, per 7T° 
trigger detected at midrapidity, in p+p, peripheral d+Au, 
and central d+Au collisions for varying trigger 7r° pt- 
Figure shows the same correlations for trigger clusters 
where the cluster-7r° pairs are both detected at forward 
rapidity. The constant pedestal, bo, was subtracted from 
the correlation function. The correlations were corrected 
for the forward ir° detection efficiency and for the combi- 
natoric background beneath the 7T° peaks in the photon- 
pair invariant mass spectra. This background is deter- 
mined by measurement of the azimuthal correlations for 
photon-pair mass selections adjacent to the 7r° mass win- 
dow and from studies of simulated jet events from pythia 
events processed through PISA. 

For the midrapidity/forward-rapidity correlations 
(Fig. [J), due to the large pseudorapidity gap of A7y~3.3 
between the hadrons, only an away-side peak (A(f>=ir) 
is seen. For the forward-forward correlations a near-side 
peak (A0=O) is also present (see Fig. [5]). The yields 
and widths of the correlated pairs are extracted by fits 
to an away-side Gaussian signal shape plus a constant 
background (bo)- The fit to the forward- forward corre- 
lations has an additional Gaussian signal for the near- 
side peak. The pedestal is determined from a fit in the 
midrapidity/forward-rapidity correlations and is consis- 
tent with the pedestal level found based on the assump- 
tion that the signal yield is at the minimum of the cor- 
relation function - zero yield at minimum (ZYAM) [20j |. 
In the forward-forward correlations the ZYAM pedestal 
is used in the yield extraction. Additional systematic 
uncertainties of up to 30% (not shown in Fig. [2]) are 
ascribed to the near-side peak due to corrections for 
resonance decays that contaminate the jet signal, and 
due to the acceptance loss around the trigger particle of 
Acf> x A77 « 0.5 x 0.5 rad, resulting from the minimum 
separation cut of one tower between cluster peaks in the 
MPC. The acceptance loss gives rise to the decrease ob- 
served for the near side peak. 

Figures [Hand [2] show that the away-side peak for d+Au 
central collisions is suppressed compared to p+p colli- 
sions and peripheral <i+Au collisions. This effect is large 
for the midrapidity/forward-rapidity correlations (Fig.[T]) 
and becomes even larger when both particles are required 
to be in the forward-rapidity region (Fig. [2]) . 

For the midrapidity/forward-rapidity correlations, 
within their large uncertainties the Gaussian widths of 
the away-side correlation peak remain the same be- 
tween p+p and central d+Au and the broadening pre- 
dicted in the CGC framework in Ref. 11] is not ob- 
served. For example, in d+Au central collisions, a — 
0.93±0.09 stat ±0.139 syst for = 1.25 GeV/c and trig- 
ger particle momentum 2.5 < p T < 3.0 GeV/c, while 
a = 0.97±0.07 stat ±0.08 syst for p+p collisions. For the 
forward- forward correlations, the measurement does not 
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FIG. 3: (color online). Relative yield JdA versus (N co i{) for forward-rapidity (3.0 < r\ < 3.8) 7r°'s paired with (left) midrapidity 
(\tj\ < 0.35) hadrons and 7r°'s and (right) forward-rapidity (3.0 < r\ < 3.8) cluster- 7r° pairs for the indicated combinations of 
Pt ranges. Also plotted as inverted solid triangles are the values of the forward n° RdA- Around each data point the vertical 
bars indicate statistical uncertainties and the open boxes indicate point-to-point systematic uncertainties. The gray bar at 
the left in each panel represents a global systematic scale uncertainty of 9.7%. Additional centrality dependent systematic 
uncertainties of 7.5%, 5.1%, 4.1%, and 4.8% for the peripheral to central bins, respectively, are not shown. The (TVcoii) values 
within a centrality selection are offset from their actual values for visual clarity (see text for actual {N co n) values). 



discern whether there is appreciable broadening between 
d+Au and p+p collisions, as the ZYAM pedestal deter- 
mination can bias the widths to smaller values. 

The observed suppression is quantified by studying 
the relative yield, JdA 2l| . of correlated back-to-back 
hadron pairs in d+Au collisions compared to p+p colli- 
sions scaled with (N co \{), 



JdA — IdA x R dA = 



1 



pair 
T dA 



l<JdA 



(N c , 



pair 
Tpp 



(1) 



where R dA 



(l/(iVcoii)) • Ka/^aVKp/ctpj,) is the 
usual nuclear modification factor for trigger particles t, 
and a, <r*, and cr pair are the cross sections (or normalized 
yields) for the full event selection, trigger particle event 
selection, and dihadron pair event selection. IdA is the 
ratio of conditional hadron yields, CY, for d+Au and 
p+p collisions: 



CY 



J rf(A0)[rfAVrf(A0) -b ] 
N { x e a x Arj a x Ap T 



(2) 



with the acceptance corrected dihadron correlation func- 
tion dN/d(A(j)), the number of trigger particles N f , the 
detection efficiency for the associated particle e a and the 
level of the uncorrelated pedestal in the correlation func- 
tions bo. The integral is taken over the Gaussian fit of 
the away-side peak. The JdA uncertainties include a sys- 
tematic uncertainty from the ZYAM pedestal subtrac- 
tion. In determining this uncertainty it was assumed that 



changes between d+Au and p+p in the Gaussian away- 
side width remain below a factor two. This upper limit 
is based on the small observed changes in width in the 
midrapidity /forward-rapidity correlations and the corre- 
lations studied previously with the PHENIX muon spec- 
trometers Q • The JdA is calculated from the measured 
IdA and R dA for the forward-rapidity trigger correlations 
with the new ir° RdAu = Rd.Au determined in the MPC. 
For the midrapidity trigger correlation s, p ublished values 
for R dA from the 2003 RHIC run [H, 13 were used. 

Figure [3] presents JdA versus (iV co n) for forward- 
rapidity 7r°'s paired with midrapidity hadrons and n s, 
and for 7r° 's and clusters paired at forward rapidity. The 
JdA decreases with an increasing number of binary col- 
lisions, (iV co n), or equivalently with increasing nuclear 
thickness. The suppression also increases with decreas- 
ing particle px and is significantly larger for forward- 
forward hadron pairs than for midrapidity/forward- 
rapidity pairs. The observed suppression of JdA ver- 
sus nuclear thickness, pr and r\ points to large cold 
nuclear matter effects arising at high parton densities 
in the nucleus probed b y t he deuteron, consistent with 
predictions from CGC [12|. This trend is seen more 
clearly in Fig. [4] where JdA is plotted versus = 
((pTi)e - ^ 1 ^ + (jpT2) e^^^/y/SNN for all pair selections 
in X] and pr- In the case of 2— >2 parton scattering, where 
two final state hadrons carry the full parton energy, 2=1, 
the variable a^u 8 wou ld be equal to (ia u ), which is the 
average momentum fraction of the struck parton in the 
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Aii nucleus. 



probing more forward rapidities. 
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4: (color online). JdA versus a;^ s for peripheral (60- 
and central (0-20%) d+Au collisions at y'sjviv = 200 
GeV. The statistical error bars and systematic uncertainty 
boxes are the same as in Fig. [3] Above x l ^ s > 10~ 3 , some 
data points were offset from their true xj^ to avoid overlap. 
The leftmost point in each group of three is at the correct 

frag 



Because the fragmentation hadrons on average carry a 
momentum fraction (z) < 1, will be smaller than 

(xau)- Based on previous studies by PHENIX at midra- 
pidity, the mean fragmentation (z) is expected to be be- 



tween 0.5-0.75 22]. In general the theoretical extrac- 



tion of xau from the measured px and r\ will differ from 
the leading order QCD picture of 2— >2 processes used 
above. Also, at modest p^'s the interpretation of the 
measured correlation functions as high energy 2— >2 par- 
ton scattering accessing low x may be limited by con- 
tributions from processes with small momentum transfer 
Q 2 . Future theoretical analysis will be necessary to eval- 
uate these and other contributions from different nuclear 
effects 4 10] on the observed large suppression in JdA- 
These analyses could additionally be complicated by the 
presence of hadron pairs originating from multiparton in- 
teractions [23| that might not probe gluon structure at 
low x Au . 

In summary, measurements of the inclusive ir° yield 
at forward rapidity, of the back-to-back correlated yield 
of cluster- ir° pairs in the forward-rapidity region, and of 
the correlated yield of forward-rapidity 7r°'s with midra- 
pidity 7r°'s or hadrons in p+p and d+Au collisions at 
^snn — 200 GeV were presented. The correlated yields 
of back-to-back pairs were analyzed for various kinematic 
selections in pt and rapidity. The forward-central pair 
measurements show no increase in the azimuthal angular 
correlation width within experimental uncertainties. The 
correlated yield of back-to-back pairs in d+Au collisions 
is observed to be substantially suppressed relative to p+p 
collisions with a suppression that is observed to increase 
with decreasing impact parameter selection and for pairs 
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